Frequency dependent admittance of a two-dimensional quantum wire 
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The frequency dependent conductance of a two- 
dimensional quantum wire is computed using a current con- 
serving formalism. The correction to the dc-conductance due 
to a time-dependent potential is related to the local partial 
density of states which we compute numerically. The current 
conservation is explicitly confirmed by computing the global 
density of states and comparing it with a quantity which is 
related to the electron dwell time. Our calculation clearly re- 
veals the physical meaning of the various partial density of 
states. 
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The physics associated with quantum conduction in 
various low dimensional systems has been a major focus 
in current sohd state research. Due to advances in con- 
trolled crystal growth and lithographic techniques, it is 
now possible to fabricate various two-dimensional sub- 
micron structures where accurate experimental measure- 
ments can be mad«l. On the theoretical side, the un- 
derstanding of quantum transport in these very SHiall 
systems has been advanced by the Landauer theorynl of 
one-dimensional transport. In cases where multi-probes 
are attached to the conductor, Biittiker has provided 
the useful formula for computing various conductancesB. 
The Landauer-Biittiker formalism has been applied to a 
large class of quantum transport problems in the meso- 
scopic and ballistic regimes where the external potential 
is static. On the other hand, when external potential has 
a time-dependent oscillating component, the Landauer- 
Biittiker formalism can not be directly applied. As shown 
by Biittiker and his co-workersu that a direct application 
of the original approach of Landauer-Biittiker formalism, 
where the internal electric potential distribution is not 
needed, can not yield electric current and charge con- 
servation. Hence the coherent quantum transport theory 
for mesoscopic and ballistic conductors must be extended 
when time- varying external potentials are to be included. 
Clearly the problem of computing frequency dependent 
conductance for mesoscopic and ballistic conductors is a 
very important subject from both theoretical and appli- 
cation points of view. 

ReceatUf in a series of articles, Biittiker and co- 
workera3Ll have developed a theoretical formalism for 
investigating frequency-dependent coherent quantum 
transport for mesoscopic and ballistic conductors when 
electric potential in the leads far away from the scatter- 
ing region of the conductor is time varying. This is an 
important advance in the theory of mesoscopic physics 



and warrant further development of the theory and ap- 
plication. In simple one dimensional systems such as a 
(5— function potential, a perfect quantum wire, or an ID 
quantum well, the frequency dependent admittance can 
be explicitly calculated since the scattering matrix can 
be obtained analyticallyEl. However, in truly 2D sys- 
tems the application of the formalism is less straight- 
forward. This is due partly to technical difficulties of 
computing certain quantities such as the various partial 
density of states (PDOS, see below). These quantities, 
which appear naturally in the theoretical formulation of 
the problemQ, have deep physical meaning which is not 
obviously clear and needs to be revealed. To the best 
of our knowledge, there are so far no studies and direct 
calculationda of the various partial density of states in 
2D. But these quantities must be obtained if the fre- 
quency dependent admittance is to be computed from 
first principles. The purpose of this paper is to report 
our investigations on the dynamic transport properties 
of a 2D quantum wire system. 

In particular, we consider a T-shaped quantum wire 
structure as shown in Fig. (1). The two probes extend 
to X = ±cxD while the scattering region is provided by the 
T-shaped junction as shown being bounded by the two 
dotted lines. Tjliis 2D quantum wire has been studied by 
many authorsOo for the case of static transport. While 
this is a simple 2D system, we found that much physi- 
cal intuition about various quantities can be obtained for 
the dynamic transport (see below). We have computed 
the frequency dependent admittance for this-system us- 
ing the current conserving formalism of Ref.u. The cor- 
rection to the dc-conductance due to a time-dependent 
potential is related to the local partial density of states 
which we compute numerically. The current conservation 
is explicitly confirmed by computing the global density 
of states and comparing it with a quantity which is re- 
lated to the electron dwell time. Our direct calculation 
clearly reveals the physical meaning of the partial den- 
sity of states which plays an essential role in the dynamic 
transport theory. 

For the sake of presentation we briefly review the cur- 
rent conserving dynamic transport formaliMH of Biittiker, 
Thomas, Pretre, Gasparian and ChristenQCl. Their the- 
ory proceeds in three steps. The first step is to determine 
the current and charge density response of the system to 
the time-variation of the external potential. This leads 
to an admittance matrix with elements given by, to first 
order in to. 



f^p{Lo)=glp{Q)-iue'(dN^p/dE) 



(1) 



1 



where the first term on the right hand side is the ad- 
mittance (conductance) when frequency a; = and this 
term is given hy the usual Btittiker multi-probe conduc- 
tance formulaB. The second term gives a correction due 
to the time-dependence of the external potential, and is 
determined by the global partial density of stateaJ (global 
PDOS), dNa/3/dE. Here the indices a/3 denote scattering 
from a probe labeled by f3 to that labeled by a. This ad- 
mittance matrix g^p{i^) does not conserve current, since 
'^adafji^) 7^ 0. To correct this problem, one must con- 
sider the internal potential distribution induced by the 
external perturbation. This is computed in the second 
step of the theoryu, and the result is rather simple if 
only Thomas-Fermi linear screening is included. In the 
final step, one computes the currents in the probes in- 
duced by this internal potential distribution, which leads 
to an additional term in the admittance matrix. Thus 
the total admittance is given by — g'^ + g^, with 
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In Eq. (^ the quantity dn{r,a)/dE is the injectiv- 
ity which measures the additional local charge density 
brought into the sample by the oscillating chemical po- 
tential at probe a. Without a magnetic field, the emissii^ 
ity dn{P,r)/dE equals to the injectivity. It was shownEj 
that these quantities could be computed from the elec- 
tron dwell time or the scattering Green function. Fi- 
nally, the total local density of states is dn{r)/dE = 
J2a dn{a,r)/dE. With these expressions, we obtain the 
final formula for the frequency dependent admittance to 
linear order in lo, 
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It is now straightforward to prove that current is 
conserved since the admittance matrix g^ satisfies 
'^adapi'^) ~ ^- This can be seen by realizing that 
J2a d^ap/dE = dNp/dE is the injectance which is iden- 
tical to J2a^al3- 

Let's apply this formalism to the T-shaped 2D conduc- 
tor. The incident electron comes from probe 1, it scat- 
ters at the T-junction, and then reflects back to probe 1 
or transmits to probe 2. To be concrete we have fixed 
the wire width to be W, the width and height of the 
side stub are fixed at aW and bW with a — b — I. 
The units are fixed by h'^/2m = 1 with m the effective 
mass of the electron, and lengths measured in terms of 
W. For simplicity we have focused on the first transport 
subband only, thus the incoming electron energy is re- 
stricted: {n/W)"^ < E < {2tt/WY. In this case elements 
of the scattering matrix s^^ are simple complex numbers. 
With more than one subbands they become matrices in 



the subband space. However the computation proceeds 
in similar fashion with one or more than one subbands. 
We have solved the quan|tum scattering problem using 
a mode matching methodllj, where the wavefunction in 
the scattering region was expanded using a suitable ba- 
sis set and in particular 50 modes were included. We 
have checked the convergence that including more modes 
in the scattering region essentially does not change the 
results. The scattering probabilities were obtained by 
matching the wavefunctions and their derivatives at the 
boundaries between the scattering region and the probes. 

To compute the admittance we need to know various 
partial density of states. These quantities turn out to 
have very interesting and physically clear property once 
they are plotted. The global PDOS is relatednk) the 
scattering matrix, and is approximately given by 
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For a finite scattering volume there are corrections to 
this expression of the order O(^) w^ere A is the electron 



wavelength and L the system sizeQ. For large system 
sizes or large energies, Eq. (|^) is adequate. Later when 
we check the current conservation, the correction term 
will be added. Since in 2D the expression for the scatter- 
ing matrix can not be written down analytically, we have 
decided to perform the energy derivatives of (|5|) numeri- 
cally. In particular we computed Sa/3{E) at various values 
of energy E and used a 5-point numerical derivative to 
find dSap/dE. The behavior of the global PDOS is plot- 
ted in Fig. (2) together with the transmission coefficient 
T2i{E). As our quantum wire system is very transmis- 
sive, T21 has large values in general except at the two 
resonance levels where T = 0. In a previous work we 
have shown that the energies where T — corresponds 
to the quasi-bound states located in the T-junctionEJ. 
The behavior of the global partial density of states coin- 
cides well with that of the transmission coefhcients. In 
particular dNu/dE, which is the global PDOS for reflec- 
tion, peaks at the energies where transmission coefficient 
is minimum or reflection maximum. On the other hand, 
dN2i/dE, which is the global PDOS for transmission, 
takes minimum value when T is minimum. Indeed, a 
general discussion for an one-dimensional system with a 
symmetric scattererlZl reveals that dNu/dE RdN/dE 
while dN2i/dE - TdN/dE with R and T the reflec- 
tion and transmission coefficients, and dN /dE the global 
DOS. Hence these quantities have vivid physical mean- 
ing. 

The next quantity of interest is Dap given by Eq. 
To compute this quantity, we first find thej-injectivity 
which is given by the scattering wavefunctionu in the T- 
junction, at zero temperature this is given by 
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where J is the incoming particle flux. Clearly a spa- 
tial ptegration of this quantity gives the electron dwell 
timellilM Tq, . In general Dap is obtained with two calcu- 
lations of the wavcfunction for particles coming from left 
and from right. Fig. (3) shows this quantity as a func- 
tion of the electron energy. Both Dii{E) and D2i{E) 
peak at the energies where the transmission takes the 
minimum values. This is understandable since tbese spe- 
cial energies correspond to the scattering statest3 where 
the electron dwell time takes maximum value, and the 
injectivity represents essentially the electron dwelling in 
an area dr at position r irrespective where it is finally 
scattered. 

With the global partial density states and the quan- 
tity Dap calculated, we have thus obtained the admit- 
tance gipiuj) from Eqs. (|l]), (||), and (||). Fig. (4) shows 
the w-dependent part of , g^p = dNap/dE ~ Dap as 
a function of energy E. Apart from the prefactor we^, 
this quantity g is the imaginary part of the admittance, 
and is called emittance. At the quasi-bound state lev- 
els where transmission coefficient T2i{E) takes minimum 
values, gap{E) also takes extremal values. It is interest- 
ing to observe that the transmissive part ^21 takes mini- 
mum value where T{E) = 0, while the reflective part gn 
takes maximum value at the same energy. Thus in this 
sense, the dynamic part of the admittance gap has the 
same behavior as the static admittance gapi^ = 0) as 
a function of energy. Furthermore, gap changes sign as 
energy is varied: the system responds either capacitively 
when gii = —§21 = —gi2 > 0, or inductively otherwise. 
Hence at the resonance (where T « for our system) §21 
and gii are capacitive. Fig. (4) shows the clear crossover 
between the capcitive and inductive responses for this 2D 
quantum wire as energy is changed. 

We can now explicitly confirm the current conserva- 
tion by summing up the admittance matrix elements and 
check whether or not J2a 9ap = 0- Since the global 
PDOS obtained using Eq. is not exact, a correctioa 
should be added. For an ID system Gasparian et. alM 
have shown that 
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where E is the electron energy. Unfortunately for a 2D 
system such as ours, the Green's function can not be 
written down analytically hence how to derive a similar 
correction term as that in (|^) is unclear. However since 
we have numerical results of all the quantities, a reason- 
able correction term can easily be obtained. We found 
that the same form as led to almost perfect current 
conservation, provided we use the transport energy fc^ as 
the energy £' in (|^). Our data in Fig. (5) clearly and 
unambiguously shows that the following is established 



E 



dN. 



With this result, we have thus explicitly shown0 the cur- 
rent conservation 17^^ = 0. Since the correction term 
is inversely proportional to fc^, it plays a role only at 
small k. This can already be seen in Fig. (4) where the 
two curves add up to zero except at low energies. 

In summary we have, for the first time, implemented 
the current conserving dynamic conductance formalism 
for 2D metallic conductors. In the simple case of ap- 
plying Thomas-Fermi linear screening for the interacting 
electrons (hence the analysis is more suitable for metallic 
samples), the current response to the internal potential 
can be computed from the electron dwell time using the 
scattering wavefunctions. The frequency dependent ad- 
mittance of a T-shaped 2D quantum wire is calculated to 
linear order in frequency. We have explicitly confirmed 
the electric current conservation and found that a cor- 
rection to the total density of states is needed at low 
energies, in the same fashion as that of ID case. It is 
very interesting to clarify the physical meaning of the 
partial DOS. These quantities provide information on 
the density of carriers and the transmission of carriers 
from one contact to another. At quantum resonances 
these quantities take extremal values. There are many 
further applications of this important theoretical formal- 
ism, to situations involving multi-mode, magnetic fields, 
and non-linear screening. Another very challenging exQ 
tension is to push the theory to higher order in frequency^ 
and formulate the theory in a way that permits a compu- 
tational implementation. We hope to be able to report 
these results in the future. 
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FIG. 1. Schematic plot of the T-shaped quantum wire. The 
wire width is W, the side- stub width and height is aW and 
bW . The two dotted lines separate the scattering region from 
the two probes. 



FIG. 2. Global partial density of states and the transmis- 
sion coefficient as functions of electron energy E. Solid line: 
transmission coefficient T21; dotted line: dNn/dE; dashed 
line: dN2i/dE. Unit of energy is h'^/{2mW^). 



FIG. 3. The current response to the internal potential, 
Dai3, as a function of energy E. Solid line: Dn; dotted line: 
D21. 



FIG. 4. The imaginary part (dynamic part) of the admit- 
tance, ga/s = dNap/dE — Da/s as a function of energy. 



FIG. 5. Explicit confirmation of the current conservation, 
Eq. (|). Solid line: dNp/dE given by Eq. (|) and /3 = 1; 
dotted line: rp/h — Dq^ with (3 = 1. 
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Fig. (1) 




-0.1 ' ' ^ ' ^ ' ^ ' ^ ' ^ ' ' 

10 15 20 25 30 35 40 

energy 

Fig. (2) 



TempleGraph 2.4, Origin: /morgoth/users/guo@celebom - Mon Jun 10 12:39:33 1996 (data file was gpdos_dat) 




TempleGraph 2.4, Origin: /morgoth/users/guo@celebom - Mon Jun 10 12:47:27 1996 (data file was dwell_dat) 




TempleGraph 2.4, Origin: /morgoth/users/guo@celebom - Tue Jun 11 09:37:47 1996 (data file was adimttance_dat) 
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